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Algorithms  for  Rsportlng  and  Counting 
Gaomstric  Zntarsactions 


by 

Jon  L.  Bant lay  and  Th.  Ottmann 


Abstract 

An  Intaraating  class  of  "Gsonstrlc  Zntsrssctlon  Problsas"  calls  for 
dealing  with  ths  pairwiss  intsrssctions  aaong  a ast  of  N objects  in 
the  plans.  These  problsas  arise  in  many  applications  such  as  printed 
circuit  design*  architectural  data  bases,  and  computer  graphics. 
Shaaos  and  Hoey  have  described  a nusaber  of  algorithms  for  detecting 
any  two  objects  in  a planar  set  intersect.  Xn  this  paper  we 
entead  their  work  by  giving  algorithms  which  count  the  nuaber  of 
such  intersections  and  algorithms  which  report  all  such  intersections. 

Mywordst  Cooputational  geometry,  geometric  intersection  problsas. 
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1.  Introduction  s ^ ” i ^ 

Many  fascinating  aspaots  of  "Gaomatric  Intaraactlon  Problaas" 
hava  baan  brought  to  light  in  tha  racant  study  of  Shanos 
and  Hoay  [1976].  Thay  invastigatad  many  diffarant  problems  • 
defined  on  sets  of  planar  objects  such  as  "do  any  tvo 
objects  intersect?". 


Thay  pointed  out  that  such  problems  arise  in  printed  circuit 
design  (do  any  conductors  cross?) , architectural  data  bases 
(are  t%iro  items  in  one  spot?) , and  operations  research 
(linear  programming  can  be  reduced  to  an  intersection  problem) . 
Shamos  and  Hoay  have  given  many  optimal  algorithms  in  their 
paper  both  for  detecting  and  for  forming  intersections  of 
many  different  classes  of  objects . 

In  this  paper  %ra  answer:  soma  of  the  open  questions  raised 
by  Shasios  and  Hoay  by  solving  problems  of  the  form  "report 
all  intersecting  pairs  of  objects"  and  "how  many  pairs 
l»^«f*«ct7".  por  example#  we  will  give  a fast  algorithm 
for  reporting  all  Intersecting  pairs  among  a set  of  N line 
segsMnts  in  the  plane.  This  problem  arises  in  printed 
circuit  board  design#  for  crossovers  must  be  placed  at 
all  such  intersecting  points.  In  this  application  (and 
many  others)  it  is  critical  that  all  such  pairs  be  reported. 

In  Section  2 of  this  paper  we  will  study  an  algorithm  due 
to  thamos  and  Hoey  for  determining  whether  any  pair  of  a 
set  of  line  segments  intersect#  and  then  generalise  their 
algorithm  to  report  all  intersecting  pairs.  In  that  section 
and  the  following  we  will  assume  that  the  reader  is  familiar 
with  Shamos  and  Hoey  [1976).  In  Section  3 we  Will  see  how 
to  modify  the  algorithm  of  Section  2 to  solve  many  other 
problems  calling  for  reporting  all  Intersecting  pairs  of 
planar  objects.  In  Section  4 we  return  to  a special  case 
of  planar  line  segments#  namely  yahmn  all  such  segments  are 
either  horisontal  or  vertical.  This  case  does  arise  in 
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applications,  and  our  algorithm  for  reporting  all  inter* 
secting  pairs  of  such  segments  is  faster  than  for  the 
general  case  (indeed,  it  is  optimal) . We  also  solve  the 
problem  of  counting  how  many  intersections  there  are  in 
such  a set.  We  give  directions  for  further  wor)c  and 
conclusions  in  Section  5. 

2.  Intersection  of  Line  Segments 

In  this  section  we  will  examine  the  problem  of  "given  N 
line  segmencs  in  the  plane,  report  all  intersecting  pairs.” 
We  will  investigate  this  problem  by  first  describing  an 
algorithm  due  to  Shamos  and  Hoey  [1976]  for  detecting  ^ any 
of  the  segments  intersect,  and  then  we  will  modify  that 
algorithm  to  report  all  intersecting  pairs.  Zn  this  paper 
we  will  not  carefully  describe  certain  important  points 
such  as  the  representation  of  line  segments  and  algorithms 
for  deciding  if  a point  is  above  or  below  a given  segment; 
we  assume  that  the  reader  is  familiar  with  Shamos  and 
Hoey  [1976],  where  these  details  are  discussed.  Throughout 
this  section  we  will  make  the  assxsnptions  that  no  segments 
in  the  set  we  are  to  process  are  vertical  and  that  no  three 
segments  meet  at  any  ona  point  — to  confront  the  details 
for  handling  these  situations  is  cumbersome  and  not  par* 
ticularly  illuminating. 

We  will  now  briefly  review  Shamos  and  Hoey's  algorithm  for 
determining  if  any  two  of  a set  of  N line  segments  in  the 
plane  Intersect.  The  basic  process  of  their  algorithm 
"draws  vertical  lines"  through  the  endpoints  of  segments 
in  the  set.  They  make  the  crucial  observation  that  the 
positions  at  which  the  different  segments  intersect  a given 
vertical  line  define  a total  ordering  on  those  segments 
(the  "above* below"  ordering) , and  if  the  segment  set  is  free 
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of  Intersections  then  the  relative  ordering  of  any  par- 
tlcular  pair  of  segments  will  be  the  same  at  all  vertical 
lines.  (Note  that  if  a line  segment  A Is  above  segment  B 
at  one  vertical  line  and  B Is  above  A at  another,  then 
they  m\ist  have  crossed,  or  Intersected,  somewhere  between 
the  two  vertical  lines.) 

Once  we  have  observed  that  there  Is  a natural  order 
relation  on  sets  of  line  segments  with  respect  to  any 
given  vertical  line  it  Is  easy  to  describe  Shamos  and 
Hoey's  algorithm.  The  main  loop  of  the  algorithm  "sweeps" 
a vertical  line  lef t-to-right  through  the  set  of  segments, 
stopping  at  each  endpoint  (this  is  implemented  algorithmi- 
cally by  sorting  the  2N  endpoints  in  an  array  and  then 
sequentially  scanning  through  it). At  each  point  during 
this  sweep  we  maintain  the  segments  which  intersect  the 
vertical  line  defined  by  the  current  x-value,  stored 
in  the  order  relation  of  the  segments  with  respect  to  the 
vertical  line.  As  a left  endpoint  is  encountered  during 
the  sweep  we  insert  the  segment  into  the  ordering  and  as 
a right  endpoint  is  encountered  we  delete  It  from  the 
ordering.  Whenever  we  insert  a segment  into  the  ordering 
we  compare  it  against  both  of  its  "top"  and  "bottom" 
neighbors  in  the  relation  and  when  we  delete  a segment 
we  compare  the  newly  adjacent  segments  — if  a given 
segment  intersects  any  segment  then  it  Intersects  one  of 

those.  Once  such  an  Intersection  Is  found  the  algorithm  reports 
it  and  halts;  if  no  such  Intersection  Is  found  then  none 
exists  among  the  segments. 

The  correctness  of  this  algorithm  has  been  proved  by  Shamos 
and  Hoey.  They  showed  that  if  the  order  relation  R among 
line  segments  is  maintained  as  a balanced  tree,  then  the 
running  time  of  the  algorithm  is  0(N  Ig  N) . We  include  a 
pseudo- ALGOL  description  of 
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Q ^ of  all  andpolnts  of  sagmants* 

storad  In  ordar  by  x^valuas 
R ♦ //  R Is  tha  ordar  ralation  of  sagmants 

currantly  axaminad 

foraach  andpolnt  p In  Q (in  ascandlng  x-ordar) 

^ p is  tha  laft  andpolnt  -at  sagaent  s then 
Insert  s In  R^ 

check  If  s intersects  tha  segments  directly 
above  or  below  it,  and  return  that  pair.if 
it  does 

//  p is  tha  right  endpoint  of  S 
check  if  tha  sagmants  diractly  abova 
azul  below  s Intarsact,  if  so  ratxurn  that  pair; 
delate  s from  R. 

Algorithm  2.1.  Dataxmina  if  N planar  segments  intarsact. 

W*  will  now  examine  tha  more  general  problem  of  reporting 
all  intersecting  pairs,  rather  than  just  saying  whether 
or  not  there  is  at  least  one  such  pair.  This  problem  was 
ixjsed  by  Shamos  and  Hoey,  who  asked  if  there  exists  an 
algorithm  to  do  this  in  time  0(N  Ig  N -f  k),  where  k is 
the  number  of  intersecting  pairs.  They  showed  this  time 
complexity  to  be  a lower  bound  on  the  problem,  we  cannot 
answer  their  question  directly,  but  we  can  modify  their 
algorithm  to  solve  this  problem  in  time  0(N  Ig  N -f  k Ig  N) . 

The  correctness  of  Shamos  and  Hoey's  algorithm  is  due  to 
the  fact  that  if  t%«j  segments  intersect  then  at  some  point 
they  must  become  adjacent  in  the  vertical  ordering.  We 
will  now  show  how  this  fact  allows  us  to  construct  an 
algorithm  for  reporting  all  intersecting  pairs  of  segments. 

We  do  this  by  "sweeping"  a line  through  the  point  set 
before) , always  maintaining  the  correct  vertical 


ord«rlnt  in  R (••  b*foT*)»  and  than  chacklng  wh«n«v«r 
•odlfylng  R to  «•«  If  navly  adjacent  aagmanta  evar  Intar- 
aact  (aa  bafora) . Thua  tha  algorithm  wa  will  praaant  la 
aubatantlally  tha  aama  aa  Shamoa  and  Hoay'a  original, 
but  modified  to  maintain  tha  correct  total,  ordering  on 
tha  aagmanta  even  after  an  intaraactlon  la  found.  Note 
that  If  two  aagmanta  are  determined  to  Intaraact  (aoma- 
whara  to  tha  right  of  tha  currant  acan  poaltlon),  than 
they  will  be  In  tha  correct  order  up  to  tha  Intaraactlng 
point  and  at  chat  point  chay  ahould  be  "awappad"  In  tha 
order.  Having  made  thla  obaarvatlon  It  la  trivial  to 
ccdlfy  Shamoa  and  Hoay'a  algorithm.  They  updated  the 
ordar  R at  the  "critical''  tlmaa  of  entering  and  leaving 
the  llna  aapianta.  Ua  will  additionally  update  R at  tha 
"critical"  time  of  aagmant  Intaraactlon.  Our  modified 
varalon  of  Shamoa  and  Hoay'a  algorithm  la  praaantad  in 
paaudo-ALGOL  aa  iUgorlthm  2.2. 


Q 4 the  s«t  of  all  andpolnts  of  segments, 
stored  in  order  by  x-values; 

R f //  The  order  relation  among  segments 

foreach  point  p In  Q (In  ascending  x-order)  ^ 

^ p Is  the  left  endpoint  of  segment  s then 
Insert  s In  R; 

check  If  s Intersects  the  segments  imnedlately 
above  and  below  It  and  If  It  Intersects 
segment  t then  Insert  the  Intersection 
point  of  s and  t Into  Q (In  x-order) 
iJjS.  H P Che  right  endpoint  of  segment  s then 
cheek  the  pair  of  segments  directly  above 
and  below  s for  intersection  and  If  they 
meet  then  add  their  Intersection  point 
to  Q (In  x-order); 
delete  s from  R 

else  //  p la  the  Intersection  of  segment  s and  t 
report  the  pair  as  Intersecting  ; 
swap  the  positions  of  s and  t In  R 

(notice  that  they  were  and  still  are 
adjacent) ; 

check  the  upper  segment  (say  s)  for 
Intersection  with  the  segment  above 
It,  the  lower  (t)  with  segment 
below  It,  and  add  any  Intersection 
points  to  Q 

Algorithm  2.2.  Report  all  Intersections  among  N planar  segments. 


Th«  corr*ecn««s  of  Altorleha  2.2  follow*  froa  th«  fact 
that  tha  total  ordaring  R of  aagmanta  is  corractly  aaln- 
talnad  at  all  tlaas.  Tha  datalls  of  tha  proof  ara  ana- 
lofoua  to  tha  arguaant  In  Thaoraa  2 of  Shaaoa  and  Roay 
(1976].  To  iaplaoMnt  tha  algorltha  afflclantly  va  can 
stora  R as  a balancad  traa  and  Q as  a haap.  As  bafora, 
va  asauM  that  thara  ara  k intarsacting  pairs.  Tha  numbar 
of  tlaas  tha  foraach  loop  of  Algorithm  2.2  Is  axacutad 
Is  axactly  2ll4>k.  Sine*  tha  total  ordar  R can  navar  contain 
Bor*  than  N sagaants  aach  oparatlon  on  R within  tha 
foraach  loop  can  b*  parfoxaad  In  0(lg  N)  tla*.  Tha  coat 
of  aach  priority  quau*  oparatlon  on  Q Is  0(lg  (2lH‘k])  • 

0(lg  N)  (sine*  k6M^)  If  Q is  raprasantad  by  a haap.  W* 
thus  sa*  that  a cost  of  0(lg  N)  Is  Incurred  at  aach  of 
tha  0(N'44c)  Itaratlons  through  tha  loop,  so  tha  total 
nmning  tla*  of  tha  algorithm  Is  0(N  Ig  N 'f  k Ig  N). 

Mot*  that  if  k Is  vary  close  to  than  the  running  tla* 
of  our  algorltha  is  actually  graatar  than  tha  0(M^)  tla* 
of  tha  ’^laiv*"  algorltha  which  checks  all  C^)  pairs  for 
Intarsactloo . 


3.  A Ganaral  Algorltha 

In  thalr  paper  Shaaos  and  Hoey  showed  how  the  algorltha 
they  glv*  for  dataetlng  intarsactlon  among  line  segaants 
can  b*  modified  to  detect  Intarsactlon  among  sets  of 
aany  different  kinds  of  objects.  In  this  section  wa  will 
show  how  our  algorithm  for  reporting  all  Intersections 
among  line  sagaants  can  b*  aodiflad  to  report  all  intar- 
sactlons  aaong  sets  of  aany  different  kinds  of  objects. 

W*  will  axplora  this  faeat  of  our  algoritiui  by  first 
aantlonlng  general  properties  of  objects  sufficient  for 
the  corractnass  of  our  algorltha  whan  applied  to  a sat 
of  such  objects,  and  than  us*  the  ganaral  construction 
to  solve  a particular  problaa. 


Algorltha  2.2  depandad  on  three  properties  of  lin*  segments 
for  its  corractnass.  It  can  also  b*  used  to  solve  Inter- 
section problaaa  on  other  objects  as  long  as  those  objects 
display  tha  following  thra*  propartiaa. 
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PI.  A vertical  llna  through  the  object  Intersects 
it  exactly  Once. 

P2.  For  any  pair  of  objects  Intersecting  the  same 
vertical  line  It  Is  possible  to  determine  algorithmi- 
cally (at  constant  cost)  which  is  above  the  other 
at  that  line. 

P3.  Given  two  objects  It  Is  possible  to  determine 
algorithmically  If  they  Intersect,  and  If  so  to 
cooipute  their  leftmost  Intersection  point  after 
some  fixed  vertical  line. 

Property  Pi  ensures  that  an  order  relation  R will  exist  for 
any  vertical  line  and  Property  P2  ensures  that  the  relation 
can  be  computed.  Property  P3  Is  used  by  Algorithm  2.2  as 
It  adds  Intersection  points  to  Q;  the  leftmost  Intersection 
point  after  the  current  scan  position  Is  ^he  one  which 
should  be  added.  So  we  see  that  If  a class  'jf  objects  C 
has  properties  PI  through  P3,  then  we  can  report  all  Inter- 
sections among  a set  of  C's  by  modifying  Algorithm  2.2  to 
read  "C”  whenever  It  reads  "segment".  The  running  time 
of  the  modified  algorithm  Is  still  0 (N  Ig  N k Ig  N) . 

An  example  of  a class  of  objects  displaying  the  above  three 
properties  are  circular  arcs  In  the  plane,  restricted  to 
exclude  arcs  which  Include  a point  with  no  derivative.  (Note 
that  an  arc  with  such  a point  can  be  represented  by  two 
arcs  without  this  property  by  "breaking"  the  arc  into  two 
at  the  place  where  the  slope  becomes  vertical) . Such  arcs 
can  be  described  mechanically  by  giving  a circle  (center 
and  radius),  two  x-values  defining  the  endpoints  of  the  arc, 
and  one  bit  saying  whether  we  are  considering  the  upper  or 
lo%fer  part  of  the  circle  within  the  specified  x-slab. 

We  have  guaranteed  that  Property  PI  Is  satisfied  by  excluding 
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arcs  which  are  both  concave  up  and  concave  down.  It  is  a 
trivial  programming  problem  to  design  an  algorithm  showing 
that  P2  and  P3  are  both  satisfied.  This  establishes  the 
fact  that  0(N  IgN+k  Ig  N)  time  is  sufficient  for  re- 
porting all  k intersecting  pairs  eimong  a set  of  N 
circular  arcs  in  the  plane. 

One  application  of  the  above  algorithm  for  determining 
arc  intersection  is  the  problem  of  Tuclidean  fixed-radius 
near  neighbors”  .In  this  problem  we  are  given  N points 
in  the  plane  and  then  asked  to  report  all  pairs  of  points 
within  some  distance  d of  one  another  by  the  Euclidean 
metric.  Notice  that  two  points  are  within  distance  d of 
one  another  if  and  only  if  two  circles,  each  centered  at 
one  of  the  points  and  both  with  radius  d/2,  intersect. 

Thus  we  can  find  all  near  neighbors  by  considering  each 
point  in  the  set  to  be  the  center  of  a circle  of  radius 
d/2  and  then  reporting  all  intersecting  circles.  (Note 
that  we  will  have  to  "break"  each  circle  into  its  top 
and  bottom  halves,  however).  This  gives  an  0(N  Ig  N + k Ig  N) 
solution  to  the  near  neighbor  problem.  This  S2une  algorithm 
can  also  be  used  to  report  all  intersections  among  a set 
of  circles  of  varying  radii. 

4.  Intersections  of  Horizontal  and  Vertical  Line  Segments 

In  this  section  we  consider  the  special  case  of  planar 
line  segment  intersections  in  which  each  of  the  N given  line 
segments  is  either  vertical  or  horizontal.  The  problem  of 
finding  all  intersecting  pairs  in  a set  of  horizontal  and 
vertical  line  segments  arises  in  many  applications.  When 
designing  an  Integrated  circuit  conductors  are  often 
restricted  to  horizontal  and  vertical  lines;detecting  all 
crossings  of  conductors  calls  for  finding  all  Intersecting 
pairs  of  vertical  and  horizontal  line  segments.  Rectllinearly 
oriented  squares  and  rectangles  are  built  from  vertical  and 
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horizontal  llna  sagmanta.  Thus,  an  algorithm  for  finding 
Intarsactlng  pairs  of  vartlcal  and  horizontal  llna 
sagmants  can  ba  usad  to  datact  all  (proparly)  Intar- 
sactlng pairs  of  squares  and  rectangles.  Finally,  wa 
will  show  how  the  "L^  fixed  radius  near  neighbors" 
problem  (which  asks  for  all  pairs  of  N points  In  the 
plane  within  some  fixed  distance  d of  one  another)  can 
be  solved  efficiently  by  this  algorithm.  If  the  distance 
la  measured  by  the  metric. 

In  order  to  simplify  the  presentation  of  the  algorithm 
and  to  clarify  the  discussion  we  first  restrict  the  Input 
to  the  case  where  no  llna  sagmants  overlap.  All  x-valuas 
of  vartlcal  lines  and  left  and  right  endpoints  of  horizontal 
lines  are  palrwlso  distinct.  We  will  later  mention  how 
to  handle  these  cases.  In  the  problem  of  Interest  we 
are  given  N vertical  or  horizontal  line  segments  In  the 
plane.  Each  vertical  line  segment  A Is  specified  by  Its 
x-eoordlnata  x(A)  and  the  y-values  of  Its  lower  and  upper 
endpoints  bot (A)  and  top (A) . Each  horizontal  llna  segment 
is  similarly  specified  by  Its  y-coordlnate  y(B)  and  by 
the  x-values  of  Its  left  and  right  endpoints  left(B) 
and  right (B) . We  will  let  M ba  the  set  of  x-ooordlnates 
M ■ {x(A)|a  vartlcaDudeft (B)  |B  horlsontaDU 
U(rlght(B)|B  horizontal). 

Note  that  M has  at  most  2N  elements. 

The  main  loop  of  our  algorithm  sweeps  a vertical  line 
from  left  to  right  through  the  set  M.  We  use  a data  structure 
R to  store  horizontal  line  segments  currently  Intersecting 
the  vertical  lines  ordered  by  their  y-coordlnates. 

Initially  R la  empty.  Whenever  a left (respectively. right) endpoint  of  a 
horizontal  line  segment  s Is  scanned  S la  Inserted  into 
(raapactlvaly  delated  from)  the  structure  R. 

When  a vertical  segment  Is  encountered  during  the  sweep 
we  check  for  Intersection  with  horizontal  line  segments 
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in  R*  N«  dnaorlb*  this  algorithm  in  paaudo-ALQOL  aa 
JU>gorithm  4 . i 


Q tha  aat  M in  aacandlng  x-ordar  (atorad  In  auch  a way 
that  for  aach  p in  Q wa  can  racogniaa  to  which 
lina  aagmant  p balonga  and  whathar  p ia 
tha  x-valua  of  a vartical  lina  aagmant  or 
tha  laft  or  right  andpoint  of  a 
horiaontal  aagmant) . 


R ♦ 


//  Tha  ordar  ralation  among  horizontal 
lina  aagmanta  (orderad  by  y-valuaa) 


foraach 


H P 


p in  Q (in  aacanding  x-ordar)  ^ 
ia  tha  x-valua  of  tha  laft  andpoint  of  a 
horizontal  lina  aagmant  s than  inaart 
S in  R. 


alaa  ^ p ia  tha  x-valua  of  tha  right  andpoint 
of  a horizontal  lina  aagmant  S than 
dalata  S from  R. 


alaa  //  p la  x-valua  of  a vartical  lina 
aagmant  S 

datarmina  A • auccaaaor  (bot(S)«R)| 

— «-  • ••  tha  laaat  y-valua  graatar  than 

or  aqual  to  bot(S)  in  Rf 

datarmina  B ■ pradacaaaor  (top(S)/R)/ 

tha  graataat  y-valua  laaa  than  or 
aqual  to  top(S)) 

foraach  horizontal  lina  aagmant  T occurring 
in  R batwaan  A and  B do  raturn  (S«T) 
aa  an  intaraaoting  pair. 


Algorithm  4.1  Baport  all  intaraacting  paira  of  N 
planar  horizontal  and  vartical  lina 
aagmanta. 
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It  Is  easy  to  see  that  Algorithm  4 . 1 correctly  finds  all 
intersecting  pairs  of  line  segments:  V7henever  a vertical 
line  segment  S is  considered  R contains  exactly  the 
horizonal  line  segments  crossing  the  vertical  line 
X « x(s)  -p  (ordered  by  y-values) . The  algorithm  then 
reports  all  pairs  (S«T  ) where  T crosses  S.  The  data 
structure  R can  be  Implemented  as  a balanced  tree.  Hence, 
the  time  to  perform  the  operations  insert,  delete, 

successor  and  predecessor  is  0(lbg  N) . Reporting  the  inter* 
secting  pairs  can  be  done  in  time  proportional  to 
their  number  if  we  use  an  appropriate  implementation  of 
balanced  trees.  Brother  (leaf-search)  trees  of  Ottmann 
and  Six  [1976]  (see  van  Leeuven  [1976]  for  an  English 
exposition) , for  example,  can  be  used  for  this  task 
where  the  leaves  of  the  tree  are  kept  in  a doubly  linked 
list.  The  structure  Q can  be  Implemented  as  a sorted 
linear  list.  Sorting  the  elements  of  M and  storing  them 
in  Increasing  order  takes  time  0(N  logN).  Hence  the  total 
performance  time  of  the  algorithm  is  0(N  log  N k)  where 
k denotes  the  number  of  intersecting  pairs. 

To  simplify  the  presentation  of  Algorithm  4.1  we  assumed 
that  no  pair  of  line  segments  share  endpoints  or  lie  on 
the  same  line,  but  these  assumptions  can  be  removed  by 
increasing  the  "bookkeeping"  performed  by  the  algorithm. 

To  handle  the  case  that  endpoints  of  segments  might 
intersect  other  segments  we  must  be  careful  to  process 
"multiple  events"  at  a vertical  line  during  the  sweep 
in  the  order  "insert  new  horizontal  endpoints",  "check 
vertical  segments",  "delete  old  horizontal  endpoints". 

We  must  also  report  any  pair  of  vertical  (likewise 
horizontal)  segments  that  meet  — this  can  be  accomplished 
before  the  main  body  of  Algorithm  4.1  is  ever  invoked.  We 
will  sketch  the  procedure  for  detecting  overlap  among 
pairs  of  horizontal  lines;  the  case  of  vertical  lines  is 


exactly  the  same.  We  first  sort  all  horizontal  segments 
y'value  and  then  consider  only  "clusters"  of  segments 
sharing  the  same  y-value.  Within  each  cluster  we  sort  the 
endpoints  by  x>value,  and  a scan  through  the  resulting 
list  can  report  all  overlapping  pairs  In  time  proportional 
to  their  number.  Neither  of  the  special  cases  we  have 
just  sketched  alters  the  0(N  Ig  N k)  running  time  of 
Algorithm  4.1. 

We  now  briefly  describe  how  to  modify  Algorithm  4.1  for 
counting  the  number  of  Intersecting  pairs  of  N horizontal 
or  vertical  line  segments  (without  reporting  them) . We 
associate  a counter  with  each  element  A In  R which 
indicates  the  number  of  elementspreceedlng  A In  R.  When- 
ever the  sweeping  vertical  line  encounters  a vertical 
line  segment  S we  determine  A*  successor  (bot(S),R)  and 
B a predecessor (top (S) fR) } the  number  of  horizontal  line 
segments  Intersected  by  S Is  the  difference  between  the 
counter  of  B and  the  counter  of  A.  This  number  Is  added 
to  the  total  number  of  Intersecting  pairs  found  so  far. 

In  a balanced  tree  counters  can  be  updated  after  an 
Insertion  or  deletion  In  time  O(log  N) . Hence,  the 
modified  algorithm  will  report  the  total  number  of  Inter- 
secting pairs,  after  0(N  log  N)  steps  which  Is  optimal. 


Finally,  we  apply  Algorithm  4.1  to  solve  the  "L^  fixed 
radius  near  neighbors"  * problem.  In  this  problem  we  are 
given  N points  In  the  plane  and  asked  for  all  pairs  of 
points  within  some  distance  d of  one  another  by  the 
metric  (l.e.  the  maximum  coordinate  metric).  The  crucial 
observation  Is  that  two  points  are  within  distance  d of 
one  another  If  and  only  If  two  squares  of  side  length  d 
and  with  sides  parallel  to  the  coordinate  axes,  each 
centered  at  one  of  the  points.  Intersect.  Hence,  for 
finding  all  pairs  of  near  neighbors  we  surround  each 
point  with  a square  centered  at  that  point  (with  side 
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length  d and  sides  parallel  to  the  coordinate  axes) . 

Then  we  use  Algorithm  4.1  for  reporting  all  pairs  of 
intersecting  squares.  This  gives  an  0(N  log  N k) 
solution  to  the  "L^  fixed  radius  near  neighbors"  problem. 
The  same  result  was  obtained  in  a totally  different 
manner  by  Bentley,  Stanat  and  Williams  [1977]. 


5.  Conclusions 


We  will  now  briefly  summarize  the  research  described  in 
this  paper.  In  Section  2 we  showed  how  Shamos  and  Hoey's 
algorithm  for  detecting  whether  any  two  of  N planar  line 
segments  Intersect  can  be  modified  to  report  all  such 
intersecting  pairs;  the  running  time  of  the  resulting 
algorithm  was  0(N  Ig  N + k Ig  N) . In  Section  3 we  modified 
the  algorithm  of  Section  2 so  it  can  report  all  intersections 
in  planar  sets  composed  of  more  complicated  objects 
than  line  segments;  the  running  time  of  the  algorithm  was 
not  changed.  We  then  examined  a special  case  of  line 
segments  (when  each  is  either  horizontal  or  vertical)  in 
Section  4 and  showed  how  all  intersecting  pairs  in  such 
a set  could  be  reported  in  0(N  Ig  N + k)  time  and 
counted  in  0(N  Ig  M)  time.  Both  of  those  performances 
are  optimal.  The  existence  of  these  algorithms  partially 
answers  a number  of  questions  posed  by  Shamos  and  Hoey. 

These  algorithms  are  particularly  interesting  because 
they  are  among  the  first  algorithms  with  complexity 
described  not  only  as  a function  of  the  problem  input 
size,  but  jointly  as  a function  of  problem  input  and 


output  sizes. 


A number  of  geometric  intersection  problems  remain  unsolved. 
The  most  outstanding  open  problem  we  have  raised  in  this 
paper  is  whether  all  k intersections  among  N line  segments 
can  be  reported  in  0(N  Ig  N 4 k)  time.  (Perhaps  one 
reason  that  our  algorithm  fails  to  meet  that  bound  is 
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that  it  raports  tha  Intersactlon  points  sortad  by  x-valua: 
this  alght  In  Itsalf  raqulra  0(k  Ig  k)  tlma.)  Tha  reader 
Intaraatad  In  further  open  problems  may  consult  the  list 
given  by  Shamos  and  Hoay  [1976].  Perhaps  some  of  the 
methods  we  have  used  In  this  paper  can  be  applied  to 
those  problems. 
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